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Abstract 

We compare decoherence induced in a simple quantum system (qubit) for 
two different initial states of the environment: canonical (fixed temperature) 
and microcanonical (fixed energy), for the general case of a fully interacting 
oscillator environment. We find that even a relatively compact oscillator bath 
(with the effective number of degrees of freedom of order 10), initially in a 
microcanonical state, will typically cause decoherence almost indistinguish- 
able from that by a macroscopic, thermal environment, except possibly at 
singularities of the environment's specific heat (critical points). In the latter 
case, the precise magnitude of the difference between the canonical and micro- 
canonical results depends on the critical behavior of the dissipative coefficient, 
characterizing the interaction of the qubit with the environment. 
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I. INTRODUCTION 



Recent years have seen significant experimental advances in manipulation of quantum 
states in a variety of physical systems 0-0] • In addition to the intrinsic interest that these 
experiments have with regard to the fundamentals of quantum mechanics, they suggest that 
a high degree of control and coherence in simple quantum systems can be achieved, perhaps 
eventually sufficient to implement a useful quantum computation || in an assembly of such 
individual units (qubits). 

Unlike a classical computer, in which the only source of errors is uncontrolled transitions 
between the states, a quantum computation is sensitive also to random changes in phases of 
the basis states. Such changes occur due to interaction of the qubit with the environment. 
They are referred to as decoherence, and the time scale over which the phase will drift by 
an amount of order one is referred to as decoherence time t^. It is advantageous to make td 
or, more precisely, the ratio t d /t s , where t s is the switching time, as large as possible. 

In many cases, the basic operations on qubits (quantum gates) can be approximated 
as evolution of certain two-level systems under an external influence (a pulse of voltage, 
current, etc.). In this paper, we concentrate on cases when the environment is comprised by 
interacting oscillators, described for brevity by a single real scalar field 0, and the interaction 
of the two-level system with the environment is linear, with the Hamiltonian of the formP] 

Hj(t) = - J d 3 xJ(x, t)0(x, t) . (1) 

In addition to Hj, the Hamiltonian of the environment contains the free Hamiltonian Hq 
and a Hamiltonian if; n t, the latter describing nonlinear interactions among the oscillators 
themselves. Our method is sufficiently general to include fully interacting oscillator envi- 
ronments, i.e. any reasonable form of H int . One motivation for doing so is the possibility 
to consider environments that are near phase transitions (or the vestiges of such in finite 
systems) . 

The "current" J(x, t) depends on the state of the two-level system and represents its 
switching history.0 This current, of course, does not have to be the usual electric current, 
although it may coincide with it in some specific cases of interactions. For example, in a 
persistent current qubit @, where the basis states differ by the value of the electric current, J 
can indeed be interpreted as the current density, and as a component of the electromagnetic 
field. 

An important example, where ([!]) applies but J is unrelated to electric current, is the 
swap gate based on two coupled quantum dots ||. Key features of this gate are as follows. 
There is an electronic spin 1/2 associated with each dot. These spins are coupled to each 
other through an exchange interaction, and the interaction Hamiltonian is given by the 



1 Bilinear and higher-order couplings in (p can be studied in a similar way. 

2 We use word "switching" to denote a controlled transition between states of a qubit, and not 
a switching on and off the interaction with the bath. The latter interaction in general cannot be 
switched off at will; see, however, remarks on the swap gate in the next paragraph. 
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energy of the singlet-triplet splitting. The exchange interaction can be switched on and off 
by varying the potential barrier between the dots; the latter is controlled by gate voltage v, 
which can be viewed as a sum of some average voltage v and a fluctuation Sv. Then, the 
singlet-triplet splitting energy can be written as 



The first term contributes to the Hamiltonian of the qubit, while the second term describes 
the interaction of the qubit with the environment. We see that this second term is precisely 
of the form (JT]), with —dET-s/dv playing the role of the "current", and Sv the role of 
the environment. For the case when fluctuations of v are the usual thermal (Nyquist) 
fluctuations, decoherence induced by Sv was considered previously in ref. |§. 

The splitting energy Ets is substantial only during a pulse of voltage that temporarily 
lowers the potential barrier. Thus, the duration of the pulse is the switching time of the 
gate. From (0), we see that in this particular case the duration of the pulse also determines 
the duration of the interaction between the qubit and the environment. 

Whenever ([!]) applies, the evolution of the field <ft from a known initial state is completely 
determined by the current J(x, t), i.e. the switching history. In other words, the gate in 
this case works as an antenna, producing a definite "radiation" state of 0. (Again, this 
"radiation" does not have to be an electromagnetic wave, but can be any kind of propagating 
excitation.) Decoherence can be associated with the probability to emit or absorb a nonzero 
number of quanta of (p. 

Typically, the initial state of is taken to be a thermal state, with probabilities of 
different energy levels given by the canonical distribution at some temperature T. In this 
paper, we want to deviate from this practice and consider a microcanonical initial state, in 
which the oscillators are constrained to have their total energy equal to some E^ There are 
several reasons why we think that this problem is interesting and potentially important for 
analysis of various qubit designs. 

First, in thermodynamics we are accustomed to canonical and microcanonical ensembles 
being essentially equivalent in the macroscopic limit. It is interesting to see if, and to what 
accuracy, the same applies to calculations of quantum coherence, which is an intrinsically 
time-dependent quantity. 

Second, some of the environments important for current qubit designs are in fact com- 
prised by relatively few degrees of freedom. Consider, for example, the swap gate described 
by eq. (§), and suppose that the pulse of voltage is delivered to the gate via a transmission 
line. Suppose further that the line is open at one end (where it attaches to the gate) and 
closed at the other (and the pulse is obtained, say, through inductive coupling of the line to 
some control circuit). For a line of length L, the number of modes significantly populated 
at temperature T is of order 



3 The current J will be set to zero at the initial moment, which can always be done by a time- 
independent redefinition of 4>. This makes choosing a microcanonical initial state for the bath 
equivalent to choosing it for the entire qubit+bath system. 




(2) 
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(3) 



For L = 1 m and T = 0.1 K, we obtain N c g ~ 10. At this point, it is actually not obvious 
that this N e ff, i.e. the number of populated modes, is what controls the transition to the 
thermodynamic limit. However, later in the paper we show that this is indeed the case. 

Now, although an ensemble of identical oscillators baths contained in identical exper- 
imental apparata may be well described by a thermal density matrix, in each individual 
experiment the bath initially has nearly fixed energy, with some broadening due to the ini- 
tial state's preparation. If the broadening is larger than the level spacing of the bath but 
smaller than the typical energy fluctuation in a thermal state, the microcanonical initial 
state is a better approximation than the canonical one. Because the level spacing decreases 
exponentially with the size of the bath, while relative thermal fluctuations only go as inverse 
square root of iV eff , we expect that such a situation will in fact be typical for relatively small, 
"mesoscopic" environments. In this case, one may wonder how much the microcanonical de- 
coherence, induced by the interaction of the qubit with the bath, differs from the thermal 
(canonical) result. In particular, one of the main goals of our project was to see if such a 
small environment can cause any significant decoherence at all. 

Finally, because decoherence is associated with the response of the environment to 
changes in the system, one may expect that anomalously large deviations from the ther- 
mal result will occur when fluctuations in the environment are large and the relaxation is 
slow, e.g. near a critical point. Our calculation lends some supports to this idea. 

Our main results are as follows, (i) If, for a microcanonical state of energy E, we 
formally define temperature T by the usual thermodynamic formula, then the expansion 
parameter that controls the difference between the canonical and microcanonical results for 
decoherence is 1/N e s, where by definition N e g = E/kgT. For a one-dimensional transmission 
line, this is of the same order as eq. (ii) We consider an expansion of microcanonical 
decoherence in powers of 1/N e g and find that the leading difference between the canonical 
and microcanonical results is formally of order l/A^fj. We present both a general formula 
for this correction, applicable for any nonlinear environment, and an explicit formula for 
an environment with Ohmic dissipation. In particular, for the case of a transmission line 
we find that already for N e g ~ 10 the environment causes significant decoherence, which 
is practically indistinguishable from the thermal result, (iii) The 1/N e s correction contains 
a term proportional to the derivative of the heat capacity of the environment with respect 
to the temperature, dCy/dT, which becomes singular near a critical point. Although in a 
finite system there can be no "true" critical singularity, a finite enhancement of dCy/dT 
remains. It is significant in this respect that our results apply to the general case of a 
fully nonlinear environment, rather than to a collection of harmonic oscillators, for which 
no critical phenomena are expected. A simple application of finite-size scaling shows that 
the critical singularity of dCy/dT alone cannot completely cancel the l/N ef[ suppression 
factor (although it can reduce the suppression considerably). However, the correction also 
depends on the dissipative coefficient, characterizing the interactions of the qubit with the 
environment, and it is ultimately the critical behavior of this coefficient that determines 
both the size and the sign of the correction at a critical point. 

In summary, while our results are somewhat inconclusive on the critical behavior of de- 
coherence (due to the lack of understanding of the critical behavior of the dissipative coeffi- 
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cient), we obtain a clear demonstration that away from criticality even a relatively compact, 
"mesoscopic" bath of oscillators, initially in a microcanonical state, induces decoherence of 
practically the same magnitude as a truly macroscopic, thermal environment. 

In the course of the evolution, the qubit and the bath exchange energy, as described by the 
interaction Hamiltonian (|1|). Moreover, because the qubit itself is controlled by some external 
means, even the compound qubit+bath system is not strictly isolated (except at times 
before and after the switching). We assume, however, that there is no additional, "direct" 
interaction of the bath with the outside world. This seems to us a reasonable assumption, 
since in most cases one will want to isolate the qubit and its immediate surroundings from 
the larger room-temperature environment. 

Our work, then, has some elements in common with the earlier work of Jensen and 
Shankar [HI on a strictly isolated small system. These authors have observed statistical 
behavior in a numerical solution of the Schrodinger equation for seven interacting spins. 
In particular, they have found that the distribution of probabilities for one spin in their 
system closely resembles the canonical distribution expected if the full spin chain were in a 
microcanonical state. In the present context, the selected spin plays the role of the qubit 
(albeit not subject to any external control), while the remaining spins play the role of the 
bath. 

Apart from the question of perfect versus imperfect isolation of the compound 



qubit+bath system, the main difference between the work of ref. |T0[ and ours is that 



[T0| compares results for a pure initial state to those for a microcanonical ensemble, for 



a compound system, while we are interested in comparing results between canonical and 
microcanonical initial states for such a system. As we will see, that latter comparison can 
be done, for a rather general case, without a recourse to numerical integrations. Instead, 
our calculation makes use of a steepest-descent evaluation of an integral relating the canon- 
ical and microcanonical averages, the accuracy of this procedure being again controlled by 

l/iV eff . 

The paper is organized as follows. In Sect. 2 we present the definition of coherence 
as a functional of the switching history, the latter being represented by the current in 
(|I|). We discuss a suitable form of the current. Although our main results are not based 
on a perturbative expansion of coherence, we pause in Sect. 3 to describe a convenient 
way to perform such an expansion, based on the coherent-state formalism. In Sect. 4, we 
compute coherence, as defined in Sect. 2, for a thermal initial state and recover some familiar 
expressions. In Sect. 5, we construct the density matrix for a microcanonical initial state. 
In Sect. 6, we compute microcanonical decoherence. Sect. 7 is a conclusion. 

In what follows we use the system of units with h = 1 and fee = 1. 



II. DEFINITION OF COHERENCE 

If we know that at some initial time t = 0, the environment started out in a definite 
quantum state |\l/(0)), we can define coherence remaining in the qubit at arbitrary time t 
in the following way. Find the final state of the environment using the evolution operator 
Uj(t,0), where the "current" J represents the switching history of the qubit. Coherence 
equals the overlap of that final with the state that would obtain if no switching took place: 
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C(t) = (*(0)\Ut(t,0)Uj(t,0)MO)). (4) 

A decrease of the overlap with time (decoherence) is due to the divergence of the evolution 
histories of the environment corresponding to different histories of the qubit. It therefore 
reflects the measuring influence that the environment had on the qubit. 

An obvious extension of this definition to the case when the state at t = is a mixed 
state with a density matrix p(0) is 

C(t) = Tr[Uj(t,0)p(0)U^t,0)] . (5) 

It is convenient to incorporate the moments of time and t in the definition of the 
current. To save notation, we describe the environment by a single real scalar 0(x, t) with 
real- valued oscillator modes ip n : 

</>(x,t) = $> n (t)V„(x) (6) 



(generalizations are of course possible). We assume that each mode couples to some smooth 
J n (r), which is zero at r < 0, switches on at r « 0, stays on a plateau until r « t, and then 
switches off, see Fig. 1. Thus, the Lagrangian of the field is 



E 

n 



+ Jn<Pn \ + U 



(7) 



where L mt describes a self-interaction. 



Jn(X) 



FIG. 1. Profile of the current J n representing the switching history of the qubit. 

The above form of the current can describe either of the following two experimental 
setups. For the swap gate described by the Hamiltonian (0), the profile shown in Fig. 1 
represents a single gate operation: both the initial and final states of the qubit correspond 
to J n = 0. So, the switching time of the qubit is the entire time t. On the other hand, if one 
basis state of the qubit corresponds to J n = 0, and another to the plateau value J n = A n , 
then the current of Fig. 1 represents two switching operations: from the first state to the 
second and back. In this case, the switching time is the ramp time of the current, r r .0 



4 And coherence defined by (|5|) coincides in this case with what is perhaps a more familiar definition: 
the value at time t of the off-diagonal element of the qubit 's density matrix, relative to its value at 
t = 0; oj r = 2it /r r acts as a frequency cutoff. 
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For this form of the current, we can relate p(0) to the density matrix at some T, < in 
the distant past as 

p(o) = E/ J (o,r i )p(r i )E/ t (o,r f ) (8) 

(since at t < Uj and Uo coincide) and also extend t in fl5|) to some Tf in the distant future. 
In this way, we obtain coherence as a functional of J n (r): 

C[J]=Tt[Uj{T f ,T i )p(T i )Ut(T f ,T i j\ . (9) 

Further, using the environment's S'-matrix 

Sj = e iHb °° T Wj(T f , Ti)e- iH{ ^ , (10) 

where Hf ree is the Hamiltonian in the absence of self-interactions (and interaction with the 
system), we can rewrite (§) as 

C[J] = Tr [SjpiSi] , (11) 

where 

Pi = e iHb - Ti p{Ti)e- iHb - Ti . (12) 

Under the usual assumption of adiabatic switching on of the interaction in the distant past, 
Pi is independent of Tj. Thus, specifying it is a convenient way to impose initial conditions. 

Eq. (|11|) is the definition of coherence that we use in what follows. We observe that 
coherence defined in this way coincides with the generating functional of the Green functions 
corresponding to the state pj. In perturbation theory, it can be computed order by order 
with the help of the Schwinger-Keldysh diagram technique. 

If we are to have small decoherence, it is natural to assume that the currents J n are 
weak. Then, provided that the field </> does not have a nontrivial expectation value, the 
leading term in Q[J] = — lnC[J] is bilinear in J: 

Q[ J ] = -^E / dtdt'J m (t)A mn (t,t')J n (t') + 0(J 3 ) , (13) 

where A mn is the full (connected) Green function of in the state pf 

A mn {t,t') = tTi {slTl^^So^} ; (14) 

(p 1 is the field operator in the interaction representation. The real part of Q, related to the 
imaginary part of A, determines the exponential suppression of coherence due to switching 
and can be called the decoherence exponent. 

For specific calculations, we will use the following expression for the Fourier transform 
of the current: 

= ^ {e im - l) exp(-|nl/2w r ) , (15) 



where u r = 27r/r r , and A n are real constants. In the limit u r — * oo, eq. (15) becomes the 
Fourier transform of a rectangular pulse: J n {r) = A n for < r < t and zero otherwise. 



7 



III. PERTURB ATIVE EXPANSION 



Although our main results are not based on a perturbative expansion, we pause here to 
outline a convenient way to carry it out. 

As we have seen, coherence naturally acquires an exponential form. So, it is convenient 
to compute (|11|) in a representation in which the trace reduces to a saddle-point integral: 
such integrals produce exponentials automatically. A good choice is the coherent-state 
(holomorphic) representation [|ll],|l!|], which we now review. (For a scattering problem with 
a large but not macroscopic number of particles and microcanonical initial conditions, the 
coherent state representation was used in ref. JI3| .) 

Any state of the environment can be represented by an anti-analytical function ip(a*) of 
the complex variable a labeling the coherent states. Action of an arbitrary operator A is 
represented by an integral of the form 



da*da 
2m i 



* a A(b*,a)^(a*) 



where A (b*, a) is the kernel of the operator A defined by 

A(b*,a) = (b\A\a) . 

A product of two operators is represented by the convolution of their kernels: 

da*da 



(A a A 2 ) (&*,&) 



2ni 



a A 1 (b\a)A 2 (a\b) 



The S'-matrix is given by fTT|^i3 



Sj(b*,a) = (b\Sj\a) = J dfad^fVcj) e B > +B f +i I Ldt 



(16) 



(17) 



(18) 



(19) 



which contains a functional integral over the field <fi as well as ordinary integrals over the 
field's boundary values 0« in and 4>f,n- The boundary terms Bi and Bf read 



n 



_I a 2 e ^nTi _ \u n $\ n + V2uj^a n (j)i >n e- 



(20) 
(21) 



The perturbation expansion for Sj is generated in the usual way via the relation (see e.g. 
ref. 0) 



Sj = exp[i J L- mt (5/i5J)dt}S'j , 



where 



S', 



id(j)fVcj) e 



Bi+Bj+i j (L bcc +J<f>)dt 



(22) 



(23) 



and L free is the Lagrangian of free oscillators. 
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The integrals in (^) are Gaussian and can be evaluated exactly at the corresponding 
saddle points. The saddle-point equation for <f> n {t) is simply the equation of motion 

n + u?J> n = J n (t) , (24) 

while the saddle-point equations for 0, and 0/ supply the boundary conditions 

u n (pi,n + i<j>i,n = V2^e- iuJnT *a n , (25) 
c^0/,n - %n = v^e^fe; . (26) 

The solution to (|2~4]) with these boundary conditions is 

<f> n (t) = -4=Ke~ iu;nt + b* n e 1 ^} + f dt'G n {t - t')J n (t') (27) 

V J 

where G n is the free causal Green function 



G n (t — t') = -J— e - iw ^\t~t'\ _ (28) 
2u>„ 



Substituting the saddle-point solution (|27| ) into (|19"D, we obtain (cf. ref. fT2|j ) 

S'j(b*,a) = expJ2Unb* n + (t/V2^) f (ftJ„(f)[fl n e^ t + ^ t ] 

+^ | dtdt' J n (t)G n (t - t')J n {t')\ . (29) 
This can be conveniently rewritten in terms of Fourier transforms of J n and G n , defined as 

/oo 
J n (T)e tnr dr , (30) 
-co 

G n {9) = / G n (r)e in ^r . (31) 

J — oo 

We obtain 

f ? ~ i r dQ, ~ ~ 1 

^(&*,a)=exp]T a n 6; + -=[ l /:K)a n + J n K)6;] + - /— G n (fi)| J n (fi)| 2 . (32) 

This can be used in eq. ( [22|) to produce a perturbative expansion for the S'-matrix. 

IV. THERMAL DECOHERENCE 



Returning to our definition of coherence, eq. ((TTJ) , we see that in the absence of self- 
interaction we would have Sq = 1 and Sj = S'j, so that 

C[J) = Tr [S' jPl ] = J ^^e-^e- b Va*, b)S'j(b*, a) . (33) 
In particular, for a thermal initial state with inverse temperature /?, 
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Pi(a*, b) = n (l - e"^") exp (<6 n e-^«) . (34) 

n 

In this case, the integrals in fl33|) are Gaussian and can be evaluated explicitly. We obtain 
C[J] = exp(-Q[J]) with 

Q[J] = ^ ji^H! [2n B K) + l]~y ^|J„(fi)| 2 ReG„(fi)| , (35) 

where Ub{oj) = [exp (/3a;) — is the Bose distribution. Eq. (j35l) is the noninteracting limit 
of the more general eq. (|T3|). 

For an interacting environment, (|34l) is still the correct initial condition for a thermal 
state, because the interaction is assumed absent in the distant past (and the interacting 
state is obtained by an adiabatic switching on of the interaction, while maintaining the 
fixed temperature 1//3). In the limit of small J, we now use eq. fllS]), according to which 
the real and imaginary parts of Q are determined, respectively, by the anti-Hermitean and 
Hermitean parts of A mn (Q) (the Fourier transform of A mn ). The anti-Hermitean part (which 
itself is an Hermitean matrix) 

A^(n) = i [A mn (o) - A^(n)] (36) 

can be expressed through the spectral density of the environment D mn in the corresponding 
channel: at uj > 0, 

A£» = A^(-w) = ^D mn {u) coth(/?o;/2) . (37) 

D mn includes effects of the self-interaction. 

The real part of Q (the decoherence exponent) becomes 



Qb[J] = E / - r r m {^)D mn {u)J n {u) coth(/3o;/2) . (38) 
mn J o 4a; 

The imaginary part is given by 

<M = - « £ / ^(n)A^(n) j B (n) , (39) 



where the Hermitean part A' mn (Q) can be expressed through the spectral density D mn via 
a dispersion relation: 



A' mn (n) = p / 

Jo 



00 gL> 
2uj 



D m n{u) , D nm (u>) 



u — ri uj + 



(40) 



P denotes the principal value. 

For a current of the form flT5D, we can introduce also another kind of spectral density, 
which takes into account the interaction of </> with the current: 

= 7T £ A m A n D mn {u) . (41) 
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We can now rewrite the decoherence exponent as 



1 r°° duo i 
Q R {t) = - / -^F(u)[l - coaut]e- u/Ur coth(/W2) . 

71 JO OJ 



(42) 



This is, of course, a familiar expression for thermal decoherence, although it is usually 
discussed for an environment comprised by harmonic oscillators. Here, we obtain it for the 
fully nonlinear case. 

For Ohmic dissipation, when 



F(oj) 



TjUJ 



(43) 



the integral in (|42| ) coincides with an integral computed by Chakravarty and Leggett 
so we can use their result to obtain an explicit functional form of Qr{€)'. 



T41 



Q R (t) = ± ln(l + u 2 r t 2 ) + 1 1n 

Z7T 7T 



(3 Tit 

— sinn — 

Tit (3 



(44) 



One should keep in mind, though, that despite this formal similarity, the macroscopic quan- 
tum coherence (MQC) problem, considered in ref. is different from ours. In the MQC 



case, transitions between basis states occur spontaneously, while in a quantum gate they are 
externally induced. In particular, the cutoff frequency u r in our case in general depends on 
the switching method. We also reiterate that in our treatment, the Ohmic form fl4*3| ) refers 
to a fully interacting environment, rather than to a collection of harmonic oscillators. So, 
for example, the dissipative coefficient rj can now depend on temperature. 



V. MICROCANONICAL DENSITY MATRIX 

In the operator language, the microcanonical density matrix for energy E can be written 



as (cf. ref. [|13 



p = M-H{H -E)=^— j d£exp[iZ(H - E)] , (45) 

in Jc 

where H is the Hamiltonian of 0, and M is a normalization factor. The contour C runs just 
above the real axis (as shown in Fig. 2 by a dashed line): a small positive imaginary part 
of £ regulates the contribution of states with large eigenvalues of H. 
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C 



Re^ 



FIG. 2. The original contour C (dashed line), used in the definition (^) of the microcanonical 
density matrix, and the deformed contour (solid line), passing through saddle points (circles). 
Locations of the four saddle points correspond to an environment with an acoustic dispersion law 



in three spatial dimensions, see eq. (53). 



In this section, we calculate A/", starting from the normalization condition 

Af- 1 



Tip = ^— f d£e-^H£) = 1 . 
2tc Jc 



(46) 



Here Z(—i£) = Tr exp(i!;H) is the thermal partition sum analytically continued to a complex 
inverse temperature — z£. 

The integrand in ([46]) has a saddle point at £ = i/3, where [3 is determined by 

E = ~\nZ{ff). (47) 

Hence, (3 is the inverse temperature related to energy E in the usual thermodynamic fashion. 
The integration contour can be deformed to pass through the saddle point. Note that it 
is essential that the contour was originally defined to run above the real axis, as the point 
£ = is typically an essential singularity of the integrand in fl46|) . Calculation by steepest 
descent in the vicinity of £ = i(3 gives 

Af = (3(2n/C v ) 1/2 exp[E(3 + In Z{(3)} , (48) 

where Cy is the field's heat capacity:0 



5 Note that using steepest descent near £ = if] implies Cy > 0, the usual condition of thermody- 
namic stability. It is curious that we have in effect derived this condition without ever referring 
directly to the second law of thermodynamics, the usual source of such inequalities. 
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C v (P) = (3 2 ^\nZ((3). (49) 

For ( f4"gp to be a good approximation, two conditions must be satisfied. First, to use the 
steepest descent, we must have 

E0 > 1 , (50) 

and, second, no other saddle point should give a contribution larger than (f48|). The left-hand 
side of (^) is our definition of the effective number of degrees of freedom, N e ff. (Recall for 
comparison that, for a collection of noninteracting classical oscillators at temperature T, 
E/T is precisely the number of oscillators). So, (|50"D is the condition that the environment 
is relatively macroscopic. 

As for the role of other saddle points, it has, strictly speaking, to be checked case by case, 
i.e. for each specific model of the environment. As an illustration, we include here results 
for two simple cases: noninteracting (linear) environments with an acoustic dispersion law, 

uj n = v s k n , (51) 

in three and one spatial dimensions. The first case can correspond for example to phonons, 
while the second to electromagnetic waves in a one-dimensional transmission line (then, 
v s ~ c). 

For linear environments, 

In Z{-ig) = - £)ln(l - e^) . (52) 

n 

So, in the first case 

lnZ(-i£) = -const x l — + 0{L 2 ) , (53) 

where the volume of the tree-dimensional region is denoted by V, and its characteristic linear 
size by L. In the second case, 

In Z(-i£) = const x y + 0(ln L) , (54) 



where L is the length of the one-dimensional region. In both (|53|) and (|5|), the constants 
are positive. Note that in the case of a transmission line, the field is the "prepotential" , 
related to fluctuations of voltage along the line, 5v(x, t), via 

Sv = Tc a ' (55) 

where C is the line capacitance per unit length. So, the correlator of 5v at coincident x will 
be Ohmic, with the "dissipative coefficient" proportional to the (C/C) 1 ^ 2 impedance of the 
line. 

We assume that in both cases the environment is relatively macroscopic, so that the 
finite-size corrections indicated in ([53]), ([54]) are negligible. Then, in the case of eq. (|53l). 
the integrand of (fP]) has four saddle-points — at £ = @, —/3, i(3, and — i(3. The integration 
contour can be deformed to pass through the first three of these, as shown in Fig. 2, and 
we find that under condition ( |50"D the main contribution to the integral indeed comes from 
the vicinity of £ = i(3. In the case of eq. (^4|) , there are only two saddle points, at £ = ±i/3, 
and only the upper one contributes to the integral after deformation of the contour. 
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VI. MICROCANONICAL DECOHERENCE 



Using the definition (|TT|) with given by eq. (|45), we obtain microcanonical decoherence 
in the form 

C[J] =M- 1 I ^ exp {-i^E + In Z(-i^-Q[t J}} , (56) 

where Q[£, J] is the thermal decoherence analytically continued to a complex inverse tem- 
perature equal to — i£. In the limit of small J n (weak decoherence), the real and imaginary 
parts of the thermal Q are given by fl38|) , fl39|) . 

Let us compare the magnitudes of different terms in the exponent of (^) on the saddle 
point £ = i/3 with (3 determined from (57). The first two terms are macroscopically enhanced: 
they are proportional to the effective number of degrees of freedom 

N cS = E(3 . (57) 

The third term, —Q[£,J], although a sum over n,m, in most cases does not have any 
macroscopic enhancement, because the couplings J n scale as l/\/N e g, and, while the diagonal 
entries of D mn are 0(1), most of the off-diagonal entries are 0(l/iV eff ). As a result, to the 
leading order in iV e fj, the microcanonical decoherence coincides with thermal decoherence at 
inverse temperature (3, while corrections are formally 0(1/ N e s). 

Even though Q[i/3, J] is not enhanced by N eS , in some cases (e.g. for Ohmic dissipation, 
cf. (0)) it grows with t (the time for which J n is on) and at large t can in principle become 
a large correction. However, for applications to qubits, we are interested only in cases when 
Qn\i/3, J], i.e. decoherence accumulated during time t, is much smaller than one. In these 
cases, corrections to the thermal result remain formally suppressed by 1/N e g. 

Let us calculate the first of these corrections. The exponent of eq. (|(J) can be written 

as 

/[£, J] = -iEZ + InZHO - Q[£, J] = /„(£) - Q[£, J] . (58) 

Here, /o(£) = /[£, 0] has a saddle point at £ = z/3, found in the previous section. The 
corresponding saddle point of the full /[£, J] is shifted to 

£ = i/3 + 6 , (59) 

with a small 5. Treating Q[£, J] in ( j58|) as a perturbation, we find 



(60) 



where primes denote derivatives with respect to £ (so that 5 is in general complex). Note 
that f{j(ifl) = —Cy/P 2 , where Cy oc iV eff is the heat capacity given by (f49|). Therefore, 
5 = 0(l/N eS ). 

The saddle-point calculation that led to eq. (|48| ) in the previous section is now modified 
in two ways. First, both the extra term in (|5S|) and the shift of the saddle point contribute 
to the saddle-point exponent. Second, they also modify the second derivative of /, which 
determines the preexponent. As a result, to the leading order in 1/Cy, we obtain 
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c ^HU&)h{^-^}' (61) 

where subscripts following commas are used to denote derivatives with respect to (3. The 
correction to the exponent is always negative, since /o,as — Cy/f3 2 > 0. Note, however, that 
although the corrections to both the exponent and preexponent are of the same order in 
1/Cy, the first is also 0(Q 2 ), while the second is 0(Q). Thus, in the most interesting to us 
limit of weak decoherence, the correction to the preexponent is more important. 

In fact, for an interacting environment, we are not really allowed to keep the correction 
in the exponent, since the higher-order terms in Q, due to the self-interaction, can give rise 
to corrections of the same order, cf. eq. (|T3|). We nevertheless retain this correction in fl6"T|) 
(and in (|63|) below) because of the traditional interest in linear environments, for which it 
is the main 0(Q 2 ) correction. 

As an example, let us take a look at eq. ( |61]) for the case of Ohmic dissipation. We 
specialize further to the large-t limit, t 3> (3, so we can use for Qr[i{3, J] the large-t limit of 
the thermal expression (fP|): 

Qr[0,J\ = j, (62) 
while Qi, which is not Bose-enhanced, can be neglected. We find 

where T — 1//3 is the temperature. The correction to the preexponent, which is the main 
correction in the limit 



1 < t/f3 < I/77 (64) 

is negative whenever Cy / r\ is a growing function of T. That is the case, for example, for linear 
environments with acoustic dispersion laws, such as those considered in the previous section. 
However, as we discuss in the conclusion, there are interesting cases when dCy/dT < 0, and 
it is in principle possible to have a positive correction to coherence. 

For the transmission line considered in the previous section, E cx T 2 , so that Cy = 
2N cfi oc T, while r] is T- independent. Then, the preexponent in fl6"3|) is equal to 1 — r/t/4/3iV eff , 
which should be compared to 1 — i]t/j3, the expansion of the exponent in the limit (|64]) . We 
see that already for N e g = 10, the correction to the thermal result is only 2.5%. 



VII. CONCLUSION 

Our main result is the calculation of a correction to the thermal result for decoherence, 
for a system interacting with a nonlinear environment that is initially in a microcanonical 
(rather than canonical) state. We expect this result to apply when a qubit interacts mainly 
with a relatively compact, "mesoscopic" environment, whose initial spread in energy is 
smaller that the typical size of 1 / yj N e g relative fluctuations characteristic of a canonical 
ensemble. 
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The correction is given by eq. ( pi] ) for the general case, and by eq. ( p3|) for an Ohmic 
environment. We see that the correction is in general of order 1/N e g and so is typically small 
already for iV e fj ~ 10, but it may be enhanced where dCy/dT diverges, i.e. in a proximity 
of a critical point. 

Mathematically, the correction to the exponent in (pl|), ( |63D results from the shift in the 
saddle-point value of £. According to (^J), such a shift can be interpreted as a change in the 
effective temperature of the environment, due to its interaction with the system. In view 
of its relation fl38|) to the anti-Hermitean part of the full Green function, the decoherence 
exponent Qr (when it is small) can be interpreted as the probability for the system to 
emit or absorb an excitation quantum, as a result of the current switching from J n = to 
J n = A n . Since in a thermal state the emission is more probable than the absorption, it 
is easy to imagine that the change in the effective temperature will be positive, leading to 
an increase in decoherence. (The emission probability is proportional to (% + 1), and the 
absorption probability to n B ; combined, the two make the coth(/3c<j/2) factor in (PB|).) For 
example, for Ohmic dissipation in the t 3> (3 limit, we have 

which indeed corresponds to an increase in the effective temperature. We recall, however, 
that the corresponding increase in decoherence is an 0(Q 2 ) effect, subleading in the limit of 
weak decoherence. 

The correction to the preexponent, which is the leading correction in the weak- 
decoherence limit, represents a different phenomenon, namely, a change in the typical size 
of fluctuations in the environment as it interacts with the system. The enhancement of the 
correction near a critical point reflects the presence of large fluctuations at T = T c . Given 
that, as a condition of thermodynamic stability, Cy > 0, and that it peaks at T = T c , 
we notice that dCy/dT < whenever T is sufficiently close to T c from above. From (|5~T|). 
we see that in this case the term containing dCy/dT is positive, i.e. it tends to suppress 
decoherence. 

For an environment of a finite size, the singularity of Cy at T = T c appears through 
Cy l dC v /dT scaling as some positive power of the total volume. An application of the 
standard finite-size scaling techniques ||15|| , together with hyperscaling, gives Cy oc L d+a ^ u = 



L 2 I V and 

1 dCy 



Cl dT 



oc L~ l l v , (66) 

T=T C 



where L is the linear size of the volume, d is the number of dimensions, a is the specific-heat 
exponent, and v is the correlation- length one. Since v > 0, (|66"D shows that the critical 
singularity of OCy/dT cannot completely overcome the macroscopic suppression (but can 
reduce it significantly: for comparison, away from the critical point, the left-hand side of 
(|66|) scales as L~ d ). However, according to eq. (|6"TD, the part of the correction that is 
proportional to dCy/dT is also proportional to the relevant dissipative coefficient, such as 
7] in eq. (§3|). In addition, there is a part of the correction containing the derivative of rj. 
Thus, it is the scaling of the dissipative coefficient that ultimately determines whether the 
correction to the thermal result can be large enough to be experimentally observable. 
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We note that the critical behavior of the dissipative coefficient determines also the lead- 
ing, thermal part of decoherence near a critical point. As far as we can tell, however, specific 
results regarding that behavior are not readily obtainable. In particular, these dissipative 
coefficients are distinct from the usual kinetic coefficients introduced in the dynamical the- 
ory since they involve a summation over the modes of the environment, cf. eq. ([411). 



Our results, then, can be taken to underscore the importance of a study of these quantities 
for different types of interacting environments. 
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